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1. Introduction 


A well-established connection exists between discrete urn models and the 
standard curves and surfaces used in computer-aided geometric design (CAGD) [1],[2], 
[3], [4], [5]. The Bezier and B-spline blending functions both model elementary 
stochastic processes, and many of the geometric properties of Bezier and B-spline 
curves and surfaces can be derived by studying these probabilistic models [4], [5]. 
Recently Barsky has introduced a new type of spline into CAGD called the 
beta-spline [6], [7], [8]* The purpose of this paper is to try to gain some insight 
into the properties of beta-splines by applying the techniques of urn models. 


2. Beta-Splines 


Beta-splines are generalizations of B-splines. They were developed in order 
to replace the somewhat artificial concept of parametric continuity by the more 
natural notion of geometric continuity. Briefly the idea is this: Two curves L(t) 

tQ<t<t;L, R(u) uq<u<u^ are said to meet with n^-order parametric continuity (C n ) 
if and only if 


dSt = dh, 

du^ | u=Uq dt^ | t=t^ 


k 0,1, »..,n 


Unfortunately this definition depends on more than just the geometry of the curves 
L(t),R(u); it also depends on the specific choice of their parametric 
representations. A linear change of parameter u=^v f >0 will not change the shape 
of the curve R(u), but by the chain rule 


k,k k.k ,k 

d_R _ ft d_R _ ft d_L , d__L 

, k | ^ . ki , ki , k | 

dv I v=Vq du I u=Uq dt I t=t 1 dt | t=t^ 

Thus R(^v) and L(t) do not meet with n^-order parametric continuity even 

though the curves R(pv),R(u) are geometrically identical. To rectify this anomaly, 

the concept of geometric continuity is introduced. 


Two curves L(t),R(u) are said to meet with linear n^-order geometric 
continuity (LG n ) if and only if there exists a constant p >0 such that 



u=u 


0 


k, k_ 
P d L 


dt 


t=t 


1 


k — 0,1, ...,n 
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It is easy to check that the notion of linear n^-order geometric continuity is 
invariant under linear changes of parameter. Of course, this concept is not 
invariant under non-linear changes of parameter. A more general notion of 
geometric continuity (G n ) and more general constraint equations invariant under 
non-linear changes of parameter are given in [8]. 

Splines have typically been defined in terms of parametric continuity, and the 
B-splines form a convenient basis for these parametric splines. The more general 
notion of geometric continuity requires us to search for a new set of basis 
functions suitable for these new types of splines. These basis functions are 
called beta-splines. We shall now use urn models to construct beta-splines and 
study their properties. 

3. An Urn Model for Beta-Splines 

Consider an urn initially containing w white balls and b black balls. One 
ball at a time is drawn at random from the urn, its color inspected, and then 
returned to the urn. If the ball was the jt* 1 white ball to be chosen, then 
^J(w+b) additional black balls are added to the urn; if the ball was the jth 
black ball to be chosen, then p“J(w+b) additional white balls are added to the 
urn. 


We now introduce the following notation: 


t 


_w__ 

w+b 


= probability of selecting a white 


ball on the first trial 


<rA | 3 ) 


s j(t) 


£ »<o 


l+i) + ...+pj -1 


Sj(f ,t) 


probability of selecting a white ball after selecting 
exactly j white balls in the first N trials 

probability of selecting a black ball after selecting 
exactly j white balls in the first N trials 


N N 

Bj(t) = B , t) 


probability of selecting exactly j white balls in the 
first N trials 


For each fixed (5 it can be shown that the functions Bq( t), . . . ,B^(t) 
are linearly independent polynomials of degree N and they satisfy the constraint 
equations(*) 


(*) 


,k N a k ,k_N 

d B . - P d B, 

j=i = 1 j 

dt k | t=0 dt k | t=l 


k = 0,1 , . . . ,N-1 


The functions B$(t) , . . . ,B$}(t) are the beta-spline basis functions. If 

P®1, these functions are the uniform B-spline basis functions and the urn model is 

the standard urn model for B-splines [1], [4]. 


Given a sequence of control points P=(Pq, . . . ,P M ) , we can use these 
beta-spline basis functions as blending functions to construct LG^"^ continuous 
beta-spline curves in much the same way that we use the uniform B-spline basis 
functions to define 0^“! continuous B-spline curves. Define the i* 1 * 1 curve 
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segment by setting 


N_ N 

B [(?,P](t) = >_ B N ([5,t-i)P i < t < i+1 

j=0 J J 

and define the beta-spline curve by setting 

B[^»P](t) = B i [JJ ,P] (t) i < t < i+1 0 < i < M-N 

From the constraint equations (*) it follows immediately that 


d k B. +1 [(},P] _ (5 k d k B.[(3,P] 

dt k | t=i+l dt k | t=i+l 


k = 0,1,...,N-1 


Thus B[fJ,P](t) 


is an continuous beta-spline curve. 


Without moving the control points, we can alter the shape of the beta-spline 
curve B[(J,P](t) simply by changing the scalar parameter (3. The effect of 
increasing (3 is to move the curve closer to its control polygon and to bias the 
curve towards its initial control points. Thus our p corresponds to Barsky’s bias 
parameter [7] • 


In table 1 we summarize those properties of beta-spline basis functions and 
beta-spline curves which are directly derivable from the beta-spline urn model. 
Many of these properties are new and are presented here for the first time. 


4. Conclusion 


Urn models can be used to construct beta-spline basis functions and to derive 
the basic properties of these blending functions and the corresponding beta-spline 
curves. This is only the beginning; much work remains to be done. Here we have 
dealt only with the simple notion of linear geometric continuity and with the most 
elementary beta parameter. Non-linear geometric continuity leads to additional 
beta parameters and to more complicated basis functions [8]* Whether urn models 
can give us any insight into these higher order concepts still remains to be 
investigated. 
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TABLE 1 - PROPERTIES OF BETA-SPLINE BASIS FUNCTIONS AND CURVES 


Urn 


Basis Functions 


Curve 


Probability 

==> > B N (t) = 1 

Distribution 

V 

.A 1 

rt 

<_ J 

|V 

o 

1* 

1 £ 


2. 


3. 


Symmetry between 
white, and 

black, 1 


==> Symmetry Formula 

B j(f * C > * B N-j(p _1 » 1 -t) 


Counting ==> Explicit Formulas 


==> Coordinate Free 

II 

==> Local Convex Hull Property 
==> Curve Symmetry 
BIjS.PKt) = B[f 1 ,P Reverse ](M-N + l-t) 
==> Recursion Formula 


4. Relationship 
between first 
N and first 
N+l picks 


. ==> Locality of End Points 

==> Explicit Recursion Formula 
N-j (see 4) 

8 » (t) = 

J r N+1 (()) 

„ «>*(N-l)/2 H 

B n ( t) - *- * L - ==> P. does not affect 

tt <r j ( P ) k 

^ 3 d Bi 1(3 ,P] k=0,l, . . . ,N-1 

dt k | t=i+l 

N t^ 

B ( t ) = ==> P. ... does not affect 

N N l+N 

TT <r i ( P ) k 

y =l 3 d B..[p,P] k=0,l , ... ,N-1 

dt k | t=i 

=> Recursion Formula ==> Geometric Construction 

R N+l, t , _ s N Algorithm 

B j (t) f j (t)B j (t) Let i<t<i+l 

+s j_ 1 (t)Bj_ 1 ( t) P^(t)=P i+j 0<j<N 

P^(t)=f*j~ L (t-i)P*j _1 (t) 

+ Sj” L (t-i)P^”[(t) 


Then Bip,P](t)=P 0 (t) 


5. Recursion Formula 


=> Polynomial Functions 

=> Differentiability 
Constraints 


==> Polynomial Spline 

t .n-i 


6. Linear Independence 
of first N+l 

o — moments 

7. Limiting Conditions 

If P = oo at most 1 
white ball can be 
selected. 


If (4=0 at most 1 
black ball can be 
selected 

8. Adding only balls 
of the opposite 
color 

9. Two urns each with 
two colors 


> N kiN 

d B i-i . P d _lx 

dt k - | t=0 dt k | t 


=> Linear Independence 
=> Polynomial Basis 


=> Lim B^(t) = (l-t) N j=0 

- o jj^O, 

N N = 

=> Lim B W (t) = t j=N 

P~ >0 = l-t N j=N-l 

= 0 j^N-l,N 

=> Total Positivity 


=> Independent 
Distributions 

B^(s.t) = B“(s)Bj(t) 


= 1 

==> Locally Non-Degenerate 

==> Local Subdivision 
Algorithm 

Tension/Bias 

==> Lim B[/J,P](t) is the 

p->~ ' 

polygon determined by 

1 P . . . P 

U * * M— N+l 

=> Lim B[fl,P](t) is the 

f ->0 1 

polygon determined by 
p p . 

N-l* ’ M 

=> Variation Diminishing 
Property 

=> Rectangular Tensor 
Product Surfaces 
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